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Abstract 



A canonical formalism for spherical symmetry, originally developed by Kuchaf 
to describe vacuum Schwarzschild black holes, is extended to include a spher- 
ically symmetric, massless, scalar field source. By introducing the ADM 
mass as a canonical coordinate on phase space, one finds that the super- 
Hamiltonian and supermomentum constraints for the coupled system sim- 
plify considerably. Yet, despite this simplification, it is difficult to find a 
functional time formalism for the theory. First, the configuration variable 
that played the role of time for the vacuum theory is no longer a space- 
time scalar once spherically symmetric matter is coupled to gravity. Second, 
although it is possible to perform a canonical transformation to a new set 
of variables in terms of which the super-Hamiltonian and supermomentum 
constraints can be solved, the new time variable also fails to be a spacetime 
scalar. As such, our solutions suffer from the so-called spacetime problem of 
time. A candidate for a time variable that is a spacetime scalar is presented. 
Problems with turning this variable into a canonical coordinate on phase 
space are discussed. 

PACS number (s): 0420, 0460, 9760L 

1. Introduction 

Canonical quantization is well-suited for the study of collapsing matter sys- 
tems. First, by quantizing both geometry and matter, canonical quantization 
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goes beyond the semi-classical approximation used, for example, in the stan- 
dard treatment of Hawking radiation]!]. Second, by working on arbitrary 
Cauchy hypersurfaces, one can study what happens to the canonical data 
inside a black hole as one approaches the curvature singularity. Third, by 
performing a midisuperspace reduction to spherically symmetric spacetimes, 
one obtains simpler models that, hopefully, one can then solve. Canoni- 
cal quantization is thus a promising method for investigating the formation 
and evaporation of black holes, and for studying the nature of horizons and 
singularities in quantum theory. 

As a first step toward obtaining a better understanding of spherically 
symmetric gravitational collapse, Kuchaf has given a detailed and elegant 
analysis of the canonical quantization of vacuum Schwarzschild black holes. 
He was able to cast the classical and quantum dynamics of primordial black 
holes into a simple and geometrically transparent form by turning the curva- 
ture radius R and Killing time T of the Schwarzschild solution into canonical 
coordinates on the geometrodynamical phase space: 

m,p; T(r),Pr(r); R(r),P R (r). (1) 

T and R thus become embedding variables T(r) and R(r) that specify how the 
Cauchy hypersurfaces are drawn in the spacetime. The canonical variables 
m and p also have a simple physical meaning: m is the Schwarzschild mass 
of the spacetime, and p is the difference between proper times at the right 
and left infinities. 

In terms of these canonical variables, the super-Hamiltonian and super- 
momentum constraints are equivalent to 

P T (r) = 0, P R (r) = 0. (2) 

The Hamiltonian, which is a linear combination of these constraints, weakly 
vanishes, implying that the canonical variables m and p are constants of mo- 
tion. The Dirac quantization of this theory is also particularly simple. Wave 
functions \I/ = \l/(m, t; T, R] satisfying the quantum version of the constraints 
(0) are independent of T(r) and R(r). Since the Hamiltonian of the theory 
vanishes, wave functions are also independent of the label time t. The final 
result: \I/ = $?(m). 

The next step is to extend the above analysis to include a spherically 
symmetric, massless, scalar field source. 



2 



To some extent, the geometrodynamics of a spherically symmetric, mass- 
less, scalar field coupled to gravity has already been worked out. Berger, 
Chitre, Moncrief, and Nutku (BCMN)0 addressed this problem in the early 
1970's. Subsequently, Unruh|| and Haji'cek|| carefully analyzed the BCMN 
model, especially in regard to black hole evaporation and the properties of 
apparent horizons in the canonical formalism. 

But in all of these treatments, the action for the coupled system is reduced 
to a privileged foliation specified by the vanishing of the "radial" momen- 
tum. The Cauchy hypersurfaces are no longer arbitrary; they are selected 
by the above slicing condition. For the vacuum theory, this slicing condi- 
tion amounts to working on the surfaces of constant Killing time T. These 
hypersurfaces thus cover only the static regions of the Kruskal diagram and 
fail to penetrate the horizon. Hajicek|| also chooses to foliate spherically 
symmetric spacetimes in such a way that the region interior to an apparent 
horizon is removed. 

This is not what we want to do. 

Rather, we want to be able to choose Cauchy hypersurfaces so that we 
can study what happens to the canonical data inside a black hole as we 
approach the curvature singularity. As such, we need our foliation to cover 
the whole spacetime; the hypersurfaces must be able to penetrate an apparent 
horizon. In addition, we want to know how the hypersurfaces are located 
in the spacetime. This means that we need to have embedding variables as 
canonical coordinates on phase space. Given that the super-Hamiltonian and 
supermomentum constraints can then be solved for the momenta canonically 
conjugate to these variables, the Dirac quantization of the theory would 
be described by wave functions satisfying first-order functional Schrodinger 
equations. In this way, we would avoid the difficulties associated with solving 
the second-order Wheeler-DeWitt equation. 

In other words, we desire a functional time formalism for our collapsing 
matter system. 

As mentioned earlier, Kuchaf succeeded to find a functional time for- 
malism for vacuum spherically symmetric spacetimes. The purpose of this 
paper is to present two attempts to find a functional time formalism for the 
coupled system, and to show how these attempts fail. Basically, the solu- 
tions suffer from the so-called spacetime problem oftime^. That is, the time 
variables that we introduce as canonical coordinates on phase space are not 
spacetime scalars, and hence fail to qualify as true embedding variables. 
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The plan of the rest of the paper is as follows: In section 2, we briefly de- 
scribe the canonical formalism for a spherically symmetric, massless, scalar 
field coupled to gravity. In section 3, we introduce the ADM mass as a 
canonical coordinate on phase space, thereby simplifying the constraints just 
as Kuchar did for the vacuum theory. In section 4, we define what we mean 
by the spacetime problem of time and show that the time variable T(r), 
originally introduced for the vacuum theory, is not a spacetime scalar once 
spherically symmetric matter is coupled to gravity. We also point out that, 
although one can introduce a new time variable T(r) in terms of which we 
can solve the super-Hamiltonian and supermomentum constraints, T(r) also 
fails to be a spacetime scalar. Finally, in section 5, we conclude by presenting 
a natural candidate for a time variable that is a spacetime scalar — the curva- 
ture time T of the general, spherically symmmetric, spacetime line element— 
and discuss the problems of turning this privileged spacetime coordinate into 
a canonical coordinate on phase space. 



Let (S, g) be a 3-dimensional, spherically symmetric, Riemannian space with 
coordinates x a = (r, 9, <p) adapted to the symmetry. The line element da on 
£ can be written as 



where dQ 2 := d9 2 + sin 2 6 d(f> 2 is the line element on the unit 2-sphere. Note 
that da is completely characterized by two functions A(r) and R(r) of the 
radial label r e [0, oo). The point r = is the center of spherical symmetry. 

Modulo boundary terms, the vacuum dynamics of the gravitational field 
follows from the the ADM action 



2. Canonical formalism 




(3) 
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(4) 



where 




(5) 
(6) 



:= P R R'-KP' K 
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are the gravitational super-Hamiltonian and supermomentum. The dynamics 
of a spherically symmetric, massless, scalar field propagating on this space- 
time follows from the action 

■= J dt J™ dr (vr0 - NH^ - N r Hf) (7) 

where 

:= ^(irV + iZy 2 ) (8) 

H? := Ticf)' (9) 

are the energy density and momentum density of the scalar field. The scalar 
field is coupled to gravity by adding the two actions: S G + S^. The total 
super-Hamiltonian and supermomentum are then 

H:=H G + H 4 ', H r :=H G + H?. (10) 

The details leading to all of the above results can be found in 0. 

Boundary terms and falloff conditions play an important role for vacuum 
primordial black holes. They play an equally important role for gravity 
coupled to a spherically symmetric matter source. But rather than write 
down the falloff conditions in all their detail, let us just state the main results. 
Namely, it is possible to choose falloff conditions on the canonical variables 
(0, 7r, A, Pa, R, Pji) and on the lapse and shift (N, N r ) at r = and r — > oo 
such that: (i) the total action S G + is well-defined; (ii) the t = const 
surfaces are free of conical singularities at r = 0; (iii) no boundary terms are 
needed to compensate the variation of the scalar field variables at r = and 
r — > oo; and (iv) no boundary terms other than 

- J dtN^M^t) (11) 

are needed to compensate the variation of the gravitational variables at r = 
and r — > oo. (Expression ( |TTD is equal to the boundary term at the right 
infinity for the vacuum theory. See for details.) 

For the boundary term (|TT1) written as above, the lapse function cannot 
be freely varied at r —>■ oo. If it were, we would find Moo(t) = 0, implying 
that spacetime is flat. We can remove this restriction on the variation of 
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the lapse by introducing the proper time at r — > oo as an additional 
dynamical variable. Since Noo = too, we can rewrite the boundary term as 

S as :=- j 'dtT^Moo. (12) 
The total action is then given by 

S:=S g + S* + S&e. (13) 
It is to be thought of as a functional of (0, 7T, A, Pa, R, Pr; N, N r ; r^). 



3. ADM mass as a canonical coordinate 



The total super-Hamiltonian and supermomentum 

1 
2 



H = -R~ l P R P K + -P" 2 AP A 2 



+ A^RR" - A- 2 RR'A' + -A~ l R' 2 - -A 

2 2 

+ l -A- 1 (R- V + P 2 2 ) (14) 

P r = P R R' - AP A ' + 7i(j)' (15) 

are complicated expressions of the gravitational variables (A, Pa, R, Pr)- We 
desire a canonical transformation to a new set of variables, in terms of which 
the constraints H = = H r simplify^ 

For vacuum spherically symmetric spacetimes, KuchafjfJ] found such a 
canonical transformation. He showed that the mapping (A, Pa, R, Pr) i— > 
(M, P M , R, P R ) given by 

M := \R^Pl ~ ^A~ 2 PP /2 + ^P (16) 

P M : = P-^-^Pa (17) 
R := R (18) 

Pr := P R -^R- 1 AP K -^R- 1 F- 1 AP A 

-R^A^F- 1 ( (APa)'(PP') - (AP A ) (RRJ ) (19) 



1 Alternatively, one may choose to solve the constraints H = = iJ r by imposing 
the coordinate and slicing conditions r = R, P\ = 0. This is the approach followed by 
BCMN||, Unruh[§, and Hajicek(| in their papers. We will not follow their approach 
here, since we do not want to restrict ourselves to a privileged foliation. 
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where 



(20) 



is a canonical transformation on the gravitational phase space irrespective of 
constraints or dynamics. As such, it remains a canonical transformation on 
the extended phase space that includes the scalar field variables (0, 7r). 

In terms of the new canonical variables, the expressions for the super- 
Hamiltonian and supermomentum simplify considerably: 

AH = -F- l M'R' - FP M P R + i (R~V + R 2 0' 2 ) (21) 

H r = P R R' + P M M' + 7T0' (22) 

where 

F=1-2M/R. (23) 

Notice that the left hand side of (|21~1) is the product of A ^ and H. We are 
allowed to perform such a scaling without changing the constraint H = 0. 

For vacuum spherically symmetric spacetimes, the canonical coordinate 
M(r) is the Schwarzschild mass of the spacetime. In fact, Kuchaf|2j ob- 
tained expression ( ]T6|) for M(r) by equating the ADM form of the spacetime 
line element (constructed from N, N r , and da) with the Schwarzschild line 
element 

ds * = _ (i _ ^jdT 2 + (l - ^ ) _1 dR 2 + R 2 dn 2 (24) 

for an arbitrary parametrization: T = T(t, r), R = R(t, r). It turns out that 
this reconstruction program for the mass also works for gravity coupled to 
an arbitrary, spherically symmetric, matter source. Instead of (123) , we have 



ds 2 = —G(T, R) dT 2 2M P,R) \ 1 dR 2 + R 2 dn 2 (25) 

where G(T, R) is in general different from 

F(T, R) := 1 — 2M(T, R)/R. (26) 

As shown, for example, by Synge0 and Thorne[|], M(T, R) equals the total 
ADM mass of the spacetime contained within the sphere of curvature radius 
R at the time T. Thus, the canonical coordinate M(r) has a good physical 
meaning for any spherically symmetric matter source coupled to gravity. (See 
also the papers by Guven and O Murchadha[[|.) 
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4. Spacetime problem of time 



For vacuum spherically symmetric spacetimes, the introduction of the Schw- 
arzschild mass as a canonical variable served only as an intermediate step. 
After carefully taking into account the boundary terms at the left and right 
infinities, Kuchaf|| subsequently performed a transformation that turned 
the Killing time T of the Schwarzschild solution into a canonical coordinate 
T(r) on the geometrodynamical phase space, and then solved the constraints. 
As mentioned in Sec. 1, the final result is extremely simple: 



For gravity coupled to a spherically symmetric matter source, the same 
transformation (modified slightly to account for the different topology of S) 
can be performed. Unfortunately, the final result for this case is not nearly 
as nice. First, the constraints H = = H r do not lend themselves to any 
obvious solution. Second, even if we could solve the constraints for the mo- 
menta canonically conjugate to T(r) and R(r), T(r) is no longer a spacetime 
scalar once spherically symmetric matter is coupled to gravity. Thus, this 
solution of the super-Hamiltonian and supermomentum constraints — even if 
it exists — suffers from the so-called spacetime problem of £«me||. 

Let us be more specific. Consider the transformation (t^, M, Pm) i— > 
(T, P T ) given by 



This is Kuchar's canonical transformation adapted to the topology £ = IR 3 . 
The mapping fl2"8|)-(p9"D is invertible: 



P T (r) = 0, P R (r) = 0. 



(27) 





(28) 
(29) 



M'(r) . 





(30) 

(31) 
(32) 



Tr 



OO 



-T'(r) 
T(oo). 



It also sends 



/ 

.70 



oo 



/ 

.70 



drP M {r) 5M(r) - M m St^ i-> 



dr Pj{r) 8T(r) 



(33) 
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modulo an exact differential. Thus, ((f), n, R, -Pr, T, Pj) is a canonical chart 
on the extended phase space. 

In terms of these new variables, the super-Hamiltonian and supermomen- 
tum are given by 

AH = F^PjR' + FP R V + i(R-V + R 2 0' 2 ) (34) 
H r = P R R' + PjV + 7T0' (35) 

where 

F{r) = l + ^r- f dfPj{f). (36) 
K(r) jo 

Although these expressions for AH and H r are much simpler than they were 
originally (see Eqs. (|i~4"D and (|15|)), it is still not obvious how to solve the 
constraints H = = H r . The culprits are the F _1 and F factors multiplying 
the first two terms of the scaled super-Hamiltonian (|34|). These factors are 
responsible for the nonlinear dependence of AH on Pj. We did not succeed 
to solve these equations for Pj and Pr on a general hyper surf ace .0 



2 If we choose to impose the coordinate and slicing conditions r = R, T' = (which 
are equivalent to the BCMN gauge conditions r = R, P\ = — see footnote 1), we can 
solve the constraints for the momenta Pj and Pr, and recover the BCMN-Unruh reduced 
Hamiltonian. The solution is 



where 



and 



\F(R) 



d 

dR 



R 



(1-1*1) , Pr(R) = -vr(R)0'(R) 



R" 



exp 



dR exp 



S^R) := R" 1 (R" 2 7T 2 + R 2 (// 2 ) 
The Hamiltonian for the reduced theory is simply 



dR -Pr(R) 











/ dR 


exp ^ 


/>•; 




la 









(37) 

(38) 
(39) 

(40) 



which agrees (up to a factor of i) with the BCMN-Unruh reduced Hamiltonian. This 
calculation just serves as a check on our results. As mentioned in Sec. 1, we prefer not to 
work on a privileged foliation. 
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But let us suppose, for the sake of argument, that we were able some- 
how to solve the constraints for Pj and Pr. Could we then claim that we 
found a satisfactory functional time formalism for a spherically symmetric, 
massless, scalar field coupled to gravity? The answer is "no." The reason 
is the following: A true embedding variable must be a spacetime scalar; it 
should not depend on the hypersurface from which it was constructed. If two 
hypersurfaces S and £' intersect at the same event 8 in spacetime, and if 
the canonical data on each of these hypersurfaces are related by the Einstein 
equations, then the values of the embedding variable at £ (obtained from the 
two sets of canonical data) must be equal. Otherwise, the embedding vari- 
able would assign different values to the same spacetime point. Since, as we 
shall show below, T(r) is not a spacetime scalar, it is not a true embedding 
variable. We do not have a functional time formalism for our theory, and 
this solution suffers from the spacetime problem of time0. 

The requirement that a dynamical variable be a spacetime scalar can be 
expressed in purely canonical language |J. Namely, a dynamical variable s(x) 
is a spacetime scalar if and only if: (i) the function s(x) is a spatial scalar; 
and (ii) the value s(x) is unchanged, modulo the constraints, if we evolve the 
canonical data with a smeared super-Hamiltonian whose smearing function 
vanishes at x. Condition (ii) is equivalent to 

{s(x),H(y)}&cx5(x,y) (41) 

where oc 5(x,y) is shorthand notation for terms proportional to 8{x,y). 
Thus, the Poisson bracket of a spacetime scalar with the super-Hamiltonian 
is weakly proportional to a 5-function. 

A fairly simple calculation shows that condition (|4lD is not satisfied for 
T(r). Explicitly, 

{l{r),K{f)H{f)} « F-\f) R'(f)5(r,f)+ oc 9(f - r) (42) 

where the coefficients of the terms multiplying the step function 0(f — r) 
are not weakly equal to zero unless the scalar field vanishes. Thus, T(r) 
is not a spacetime scalar, and this solution of the constraints — even if it 
exists — suffers from the spacetime problem of time. 

To conclude this section, we point out that, modulo certain technical 
difficulties^, we can perform a transformation to a new time variable T(r) in 

3 These amount to the non-invertibility of the transformation (|43|)-({46|). We can only 
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terms of which we can explicitly solve the constraints. The transformation 
(r^, M, P M , R, P R ) ^ (T, P T , R, Pr) is 

T(r) := ^ - f df F(f)P M (f) (43) 

R := R (45) 
Pr := Pr + \f P M (P" 1 - 1 + In |P|) (46) 



2' 



where 



P=1-2M/R. (47) 
The transformed constraints are 

APT = P T R + PrT' + - (r^V + R 2 /2 ) 

_I (R2 _ T '2) (>-i _ i + i n |ir|) (48) 
H r = P R R' + P T T' + 7T</>' (49) 



where 



|F(r)l = exp ( rJo /o" df PT(f) ) ' (50) 



The solution of the constraints is 



where 



Mr) = ^(^^I^Ol) (51) 
Pn(r) = -^y(PT(r)T'(r) + 7r(r)0'(r)) (52) 



\F(r)\ = R-\r) exp ( - £ «S ) £ df R'(f) exp ( £ <S^ (53) 



recover the absolute value of F if we try to invert the the transformation. (See Eq. (J5Cj) .) 
If we try to avoid this problem by restricting ourselves to one sign of F, say F(r) > 0, 
then we lose the hypersurfaces which penetrate an apparent horizon F(r) = 0. 
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and 

S^(r) := -2 R-\R' 2 - T 2 )' 1 (V tt0' - 1r' (r~ V + R 2 0' 2 )) . (54) 

If we impose the coordinate and slicing conditions r = R, T' = 0, our solu- 
tion again reproduces the BCMN-Unruh reduced Hamiltonian. (See footnote 
2.) 

Unfortunately, just like T(r), T(r) is not a spacetime scalar: 

{T(r), A(f)#(f)} « R'(f) <f(r, f)+ oc 6(f - r) (55) 

where again the coefficients of the terms multiplying the step function are 
not weakly equal to zero. Thus, this explicit solution of the constraints also 
suffers from the spacetime problem of time. 



5. Discussion 

The time variables T(r) and T(r) that we introduced as canonical coordinates 
on phase space both failed to be spacetime scalars. As such, they did not 
qualify as true embedding variables. It is important to stress, however, that 
the two attempts presented in this paper do not constitute a proof that a 
functional time formalism for spherically symmetric matter systems coupled 
to gravity does not exist. In fact, as we shall argue below, our current belief 
is that a functional time formalism for these systems does exist. We need 
only be more clever in our choice of time variable.^ 

Indeed, a natural candidate for a time variable that is a spacetime scalar 
is the curvature time of the general, spherically symmetric, spacetime line 
element 

ds 2 = —G(T, R) dT 2 + F(T, R)- 1 dR 2 + R 2 dVl 2 (56) 

where 

F{T,R) := 1-2M(T,R)/R. (57) 

(See also the discussion at the end of Sec. 3.) By its definition, the curvature 
time T is a spacetime scalar, and like the curvature radius R, T has an 
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A functional time formalism for spherically symmetric spacetimes has also been dis- 



cussed by BrahamjlO . It appears that his solution of the constraints also suffers from the 
spacetime problem of time. The "embedding" variables given in JTo[ are not spacetime 
scalars. 
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invariant geometrical meaning: (i) the surfaces of constant T are orthogonal 
to the lines of constant R, 6, and 0; and (ii) the labeling of the T = const 
surfaces is specified (up to the choice of time origin) by requiring that T 
measure proper time at R — 0. Requirement (ii) imposes the boundary 
condition G(T, R = 0) = 1 on G. 

The problem is how to turn this privileged spacetime coordinate into a 
canonical coordinate on our phase space. 

For vacuum spherically symmetric spacetimes, there is no problem. Fol- 
lowing the reconstruction program for the mass described in Sec. 3, one finds 

— T' = F- x R- x KP h (58) 

where 

/R'\ 2 /p A x 2 

(59) 




R 

As shown in 0, —T'(r) is the momentum canonically conjugate to the Schw- 
arzschild mass M(r). Then, by carefully taking into account the boundary 
terms at the left and right infinities, one can perform another transformation 
that turns T itself into a canonical coordinate on phase space. (See || for 
more details.) 

For gravity coupled to a spherically symmetric matter source, things are 
not so simple. Equation (0) is replaced by 

- G^T' = F^R^APa (60) 



where F is given by our old expression (|59|) . Thus, we have only been able to 
reconstruct the product G^T' in terms of the original gravitational variables. 
To obtain an expression for G or T' separately, we must somehow involve the 
matter variables. 

An idea that immediately suggests itself is to use one of the Einstein 
equations M: 



8tt j* dR R F~ l (T, R) (/i(T, R) + p(T, R) 



(61) 



G(T, R) = F(T, R) exp 
where 

/x:=-T T T , p:=Tg (62) 

are two components of the energy-momentum tensor T Q/ 3 for the spherically 
symmetric matter source. The problem with this approach is that the integral 
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in (|6T| ) is over a T = const surface. Even though it is possible to express 
the integrand of ( |B~lD in terms of the original canonical variables, we still 
have to evolve the canonical data from E to the T = const surface before 
we can do the integration. Since E is an arbitrary spherically symmetric 
hypersurface, E need not agree with the T = const surface anywhere. The 
resulting expression for G, and hence for T, would be non-local in time as 
well as in space. 

Another approach, which appears to be more promising, has a somewhat 
different starting point. The idea is to first reduce the Einstein-Hilbert ac- 
tion to spherically symmetric spacetime metrics of the form (p6|) and (|5T|), 
and then parametrize the resulting action to introduce the curvature time T 
and its conjugate momentum as canonical data on arbitrary, spherically sym- 
metric, hypersurfaces. In this manner, we would succeed in promoting both 
the curvature radius R and curvature time T to canonical coordinates R(r) 
and T(r) on phase space. The spacetime problem of time would thereby be 
avoided. But a possible problem with this approach is the existence of second 
class constraints. In the process of eliminating the second class constraints 
prior to quantization, we may lose T{r) as one of our canonical variables. We 
are currently investigating these issues. 
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